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Abstract 

We present a new method to study the ground state of quantum spin systems 
using the Monte Carlo techniques together with restructured intermediate states 
which we proposed previously. 

Our basic idea is to obtain coefficients in the expansion of the ground state from 
measurements of the normalized frequencies of each state on the boundary in the 
Trotter direction. The restructuring method we have proposed makes it possible to 
carry out simulations within available resources of computers. 

For a concrete example we study ladder spin systems with antiferromagnetic 
intra-chain nearest neighbor interactions, whose ground states may be equivalent to 
that of the spin 1 chain. For the ladder with ferromagnetic rungs we successfully 
show how to determine the magnitudes and relative signs of the coefficients in the 
expansion. When the ladder has antiferromagnetic rungs we examine White's claim 
that the diagonally situated pair of spins are mostly triplet states in the ground 
state. Our results strongly suggest that the system has the hidden order which is 
constructed by the specially aligned triplet states of diagonal pairs. 



Section 1 Introduction 

Recent development of experiments on the condensed matter has presented us 
very interesting systems where strong quantum effects will be realized. Among them 
are the quantum spin systems in low dimensions. Much theoretical work has been 
done to obtain qualitative and quantitative properties of these quantum effects [Q. 
One powerful tool to numerically investigate them is Monte Carlo approach using 
the Suzuki- Trotter formula [||. In a previous paper 0, where we proposed the 
restructuring method, we pointed out that proper choice of complete sets which 
describe the intermediate states in the Suzuki- Trotter formula is quite important 
in numerical calculations. Then we applied our method to a quantum spin chain 
system to show that by employing complete sets different from the conventional ones 
we can obtain improved numerical results. An essential point of our method is to 
rearrange spin states on neighboring two sites into the singlet and the triplet states 
when we construct the complete sets. 

In this paper we discuss another merit of our method which comes out in study of 
the system's ground state [^]. In order to provide a concrete example we apply the 
method to quantum spin systems on a simple ladder which have attracted many 
people for two reasons. One reason is that this model is an intermediate one between 
one spin 1 chain which has a finite energy gap (Haldane gap) and two spin 1/2 chains 
which are gapless. Another reason is that the corresponding material exists and this 
would show the high temperature superconductivity upon some doping. 

An outline of our basic idea is given in the next section together with the model 
we study as an example. In section 3 we present numerical results and the final 
section is devoted to discussions. 



Section 2 Model and Basic idea 

The Hamiltonian we study is 



H = -— ^[ffa,iffa,i+l + ^b,i^b,i+l] - IT X! ^a,i^b,i, (1) 
i=l ^ i=\ 

where Ba{b),i denotes the Pauli matrix sitting on the i-th rung of the leg a(6) and 



aa{b),N+i = (^a{b),i, N being the total number of rungs. Let us set the intra-chain 
couphng J = —1 throughout this paper and investigate the model with both positive 
and negative inter-chain couplings J'. 

The basic idea is quite simple. In the system's partition function Z its ground 
state I G) dominates when the inverse temperature (3 is sufficiently large, 

Z = trie-"'') ~ e-^^«tr(| G){G |) = e'^^^ ^(a | G){G \ a), (2) 

a 

where Eg denotes the energy of the ground state and {| a)} is a complete set. 
Accordingly we expect that the partition function Z.^ in the Suzuki- Trotter formula 
with finite Trotter number rit would be 

Zn. = EE---E(«ik"''^/"M«2)---KJe-^^/-|«i) 

^ |G)e-^^«(G|ai) (3) 

for large Ut and [3. Let us expand | G) by {| a)}, 

|G) = ^c„|a), (4) 

a 

where the coefficients c^'s should be real numbers in our representation. 
With this expansion reads 

^n.-e-/^^-Ec^.=e-^''"E4- (5) 

dl a. 

This means that c\ should be the probability to find state | a) on the boundary in 
the Trotter direction. 

Since the number of | q;)'s, which totals to 2^^ on a ladder with rungs, is 
enormously large for large it is technically important to choose a proper complete 
set for which many of Cq,'s vanish. In the conventional approach the complete sets {| 
«n)} = 1, • • • , "^i) in (0) are constructed from eigenstates of cr^({,) j with eigenvalues 
±1 denoted by -^aih),i or by Sa{h),i- The complete sets should be therefore 

\otn) =\ Sa,l, 5^,2, Sb,2, " " " Sa,Ar, Sfo^Ar). (6) 

In our approach we rearrange states Sa,i and s^^i into a singlet state | 0j) and three 

triplet states | Ij), | ®i) and | — Ij), 

. . 1. 



I ®«) ~ (I +a,i; I +a,i; ~b,i)) 5 

which we denote by Si. The complete sets we employ are 

I ttn) =1 Si,S2, - ■ ■ Sn)- (8) 

We will see later that this choice effectively reduces the number of non-zero Cq,. 

What we can measure in simulations with Z„j is not itself but c^. One may 
have a question if there is a way to know the relative sign of each Cq,. The answer is 
yes in the quantum spin case. All we have to do is to change the boundary condition 
in the Trotter direction from the periodic one to a combination of an open boundary 
I a) and a fixed boundary | ao). The corresponding partition function Z'^^ is 

Z'nt = E E ■ ■ • E(«o I e-^^/- I a,)... (a., | e"^^/- | a) 

^ ^(ao I G)e-^^«(G | a) = e'^^^ (9) 

a a 

Thus we will obtain Ca/cao = CaoCa/Ca^ in simulations with Z^^. 

With the normalization condition = 1 it would seem that the results with 
Z'^^ are enough to determine both the magnitude and the sign of the coefficients 
Cq,'s but, as we will see in the next section, simulations with are quite useful to 
obtain reliable values of magnitude with small statistical errors. 



Section 3 Numerical results 

In this section we study a ladder spin system with the Hamiltonian (|l]) using 
Monte Carlo techniques. Exact values on a small ladder are also presented to confirm 
the Monte Carlo results. 

First let us work with positive J's, which mean ferromagnetic rungs. As we 
stated in Section 2 we expect that by applying the restructuring method we can 
do with less number of non-zero coefficients in the expansion than using the 
conventional complete sets. In order to see this in a quantitative manner let us 
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examine distributions of on an = 4 and J' = 1 ladder calculated by the exact 
diagonalization. With the conventional complete sets the number of non-zero 
c^'s is sixty-eight in total, twenty of which are greater than 10~^ and sum over these 
twenty c^'s amounts to 0.808. With the rearranged complete sets (|^), on the other 
hand, total number of non-zero coefficients reduces to thirty-eight and among them 
c^'s that are greater than 10~^ are fifteen in number and 0.843 in their sum. 

In addition, for the system under consideration it is believed that the ground 
state is essentially equivalent to that of spin 1 chain, which is made of triplet states 
on sites and seems to have the hidden order that demands the | 1) and | —1) states 
should appear alternately along the chain with arbitrary number of | ©)'s between 
them ip. The complete sets (|^) are much more suitable to describe such a ground 
state. The exact results on the = 4 and J' = 1 ladder stated above, for example, 
show that all the fifteen states in (^) whose c^'s are greater than 10~^ are of this type 
while others are not. For an = 8 ladder we do not carry out exact calculations 
any more, but it is easy to see that we have only 255 candidate states in our choice, 
which can be sorted to 22 states if we take account of the translational invariance 
due to the periodic boundary condition in the space direction. Contributions of 
these particular states on larger ladders will be examined in the following Monte 
Carlo study. 

Now we present numerical results with the partition function and the complete 
sets (P). In this case we do not need to worry about the negative sign problem. We 
measure frequencies of | q;)s' appearance on the boundary in the Trotter direction 
as well as system's energy 

E = --^\nZ. (10) 

In Figure 1 we plot —E per site on an = 8 and J' = 1 ladder for several values 
of /? as a function of with Trotter numbers ranging 16 < < 32. Since Ut 

dependence of physical quantities are known to be 1/nf we make a linear fit for 
data with each value of f3 to obtain the energy in the — oo limit. We observe 
the extrapolated values almost coincide, which indicates that these values of (3 are 
large enough to realize the ground state. In Figure 2 we compare Monte Carlo 
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results on normalized frequencies measured on an = 4 ladder at (3 = 8 with 
exact numerically obtained by the exact diagonalization. We see the agreement 
is satisfactory at this value of j3. Figure 3 shows sum of contributions of the states 
with the hidden order on an = 8 ladder at /5 = 8 for several values of J' in 
the rit oo limit. We observe the sum, which indicates sum of for these state, 
increases toward unity as we change J' from 0.5 to 5. Data for values of = 4, 12 
and 16 with J' = 1 are also plotted in Figure 3 to show the dependence on A^. Since 
the results for A^ = 8 show no drastic change as a function of J' our data support 
that this system and the spin 1 chain are equivalent in properties of the ground 
state. The results also indicate that the hidden order dominates if J' is quite large, 
while it is less clear for smaller values of J'. 

Next we present results on the coefficient itself measured with the partition func- 
tion instead of (||). Note that in this case, because of the emergence of negatively 
signed weights, we should measure a physical quantity A in each thermally equili- 
brated configuration together with its sign and obtain its expectation value by 

{A) = p^ = (11) 

where Zj^iZS) is number of configurations with positive (negative) weight and 
Aj^i^AS) denotes the contribution from positively (negatively) signed configurations. 
Because of the cancellations the statistical error is larger for smaller r ratio, which 
is defined by 

Z+-Z^ 

r = — . 12 

Here we limit ourselves to an A^ = 4 and J' = 1 ladder and set the fixed boundary in 
the Trotter direction | ao) =| ©, ©, ©, ©). We then measure normalized frequencies 
of I a)s' appearance on the open boundary in the Trotter direction, /q,, following to 
(|ll|). Since /q, should be proportional to CqCq,,, we obtain the value of (?^^ from 

E(f)' = E(^)' = E(^]' = ^ = ^. (13) 

using the normalization condition I]c^ = 1. We then calculate Cq,, not c^, by 
multiplying wc^ to each fa/fao-, where we assume Cqq is positive. Figure 4 plots 
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the results together with the exact Cq,'s. We observe that the data provide relative 
signs of Ca's to Cap completely correctly in spite of rather poor value of the r ratio 
(which is around 0.15) because Z_{Z^) = for positive (negative) c^. The statistical 
errors on the magnitude of Cq, is, however, much worse than the ones measured with 
the partition function (|^). 

Finally we study the system with antiferromagnetic rungs, namely with negative 
J' in (|l]), which became much interesting recently because of White's conjecture 
0. In the limit of large | J' \ this model has the ground state that is a cluster of 
singlet at rungs as well as the excited state including a triplet state at some rung, 
with which the system has a finite mass gap. Since many studies suggest this 
mass gap survives for finite | J' |, one is tempted to think the system might be 
equivalent to the spin 1 chain. In his study based on the density renormalization 
group equation. White claimed that it is the case. It is reported that roughly 95% of 
the diagonal pairs on a ladder are triplet states in the system's ground state. Here 
we examine the model by employing the complete sets constructed with eigenstates 
of the diagonal pair. 



where S', is one of 



an) =1 S[, S2, ■ ■ ■ S'^) , (14) 



©i) — (I +a,i, —b,i+l)— I +a,i, —b,i+l)) , 



I ©i) = (I +a,i, —b,i+l)+ I +a,i, —b,i+l)) , 

\ = \ —a,i,—b,i+l)- (15) 

Using these complete sets with the partition function (^), we measure sum of nor- 
malized frequencies of the states which have the hidden order. Unfortunately we 
find that serious cancellations prevent us from obtaining statistically meaningful 
data except for small ladder size or small values of f3 shown in Figure 5. In the 
figure we see, however, that data for a = 4 ladder are enough to indicate this ratio 
increases rapidly as (3 grows and at /3 = 5 it is in good agreement with the exact 
value shown by a dashed line. Rapid increments as a function of f3 are also observed 
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for = 8, 12 and 16 ladders as far as we could measure. So we expect similar 
saturations will take place for larger ladders with large values of /3. We therefore 
conclude that the hidden order in this model is numerically observed. 

Section 4 Discussions 

In this paper we present a new method to study the ground state of quantum 
spin systems which uses the Monte Carlo techniques applied to the Suzuki- Trotter 
formula. Studying the ground state by the Monte Carlo approach will help us to 
understand its physical meanings, to make effective models and to easily measure 
some physical quantities. 

Our basic idea is to observe the states emerging on the boundary in the Trotter 
direction and count number of appearances of each state so that we can obtain 
the normalized frequencies which are related to the coefficients in the expansion of 
the ground state. This method is based on the restructuring method we proposed 
in a previous paper, which would enable us to effectively reduce the number of 
candidate states we should consider. We would like to stress here that it is essential 
from technical points of view to choose some proper complete sets, although any 
complete sets will do us from theoretical points of view. 

In order to make our discussions clear we study ladder spin systems. For the 
ladder whose rungs are ferromagnetic we show it is possible to describe its ground 
state using the Monte Carlo results obtained in our method. For the ladder with 
antiferromagnetic rungs our results support White's report that most diagonal pairs 
on the ladder are the triplet states. Our data also strongly suggest that the hidden 
order exists. 

Our method would also be applicable to other systems such as spin systems on a 
square lattice. In order to construct proper complete sets, however, one should have 
physical insights to the system under investigation. 
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Figure Captions 
Figure 1 

System's energy £" as a function of the Trotter number rif on an A^=8 and J' = 1 
ladder. The ordinate and the abscissa represent —E per site and l/nf, respectively. 
Each datum is an average of twelve measurements with different random number 
sequences, for each of which we generate 10000 configurations for the thermalizaion 
and more 10000 configurations for the measurement. Statistical error of the datum 
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is calculated from the standard deviation among these twelve runs. Here we do 
not show statistical errors explicitly because all of them are within symbols. Least 
square fitted lines for each value of (3 are also shown. 
Figure 2 

Monte Carlo results on the normalized frequencies of each state | a) in obtained 
on an = 4 and J' = 1 ladder with j3 = 5 and = 32, together with exactly calcu- 
lated non-zero c^, where coefficients in the expansion of the ground state (^). 
The abscissa represents the state | a) with non-zero coefficients which are numbered 
from 1 to 38. In each of twelve random number sequences we discarded first 2000 
configurations and used the following 10000 configurations for the measurement. 
Figure 3 

Sum of the normalized frequencies on the states possessing the hidden order in 
which is plotted versus J' for an = 8 ladder and versus for a J' = 1 ladder. Each 
datum is the extrapolated one from data for rit = 16, 20, 24 and 32 to rif oo. Data 
for finite Ut are obtained with twelve random number sequences in each of which we 
used 20000 configurations with 20000 preceding ones. Errors for the extrapolation 
calculated from the statistical errors are within symbols. 
Figure 4 

Monte Carlo results for the non-zero coefficients Cq,'s on an = 4 and J' = 1 
ladder with P = 5 and rit = 32. We calculated c^'s using the normalized frequencies 
obtained from twelve simulations in each of which 10000 configurations are observed 
after 2000 initial sweeps. The abscissa is for the state numbers which are same to 
those in Figure 2. We also plot exact values for comparison. 
Figure 5 

Sum of the normalized frequencies on the states possessing the hidden order in ([I4D 
measured on a J' = —1 ladder with the Trotter number = 32, which is plotted 
for A^ = 4, 8, 12 and 16 as a function of (3. For each (3 and A^ we carried out twelve 
simulations with different random number sequences in each of which we updated 
the system 20000 times for the thermalization and then did additional 20000 updates 
for the measurements. The r ratio ranges from ~ 0.76 to ~ 0.01. The dashed line 
is the exact value for an A^ = 4 ladder. 
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